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O Abstract 

, The first heat kernel coefficients are calculated for a dispersive ball 

whose permittivity at high frequency differs from unity by inverse powers 
of the frequency. The corresponding divergent part of the vacuum energy 
of the electromagnetic field is given and ultraviolet divergencies are seen 
to be present. Also in a model where the number of atoms is fixed the 
^ ■ pressure exhibits infinities. As a consequence, the ground-state energy for 

a dispersive dielectric ball cannot be interpreted easily. 



The ground-state energy for a dielectric ball shows ultraviolet divergencies 
still lacking physical understanding. This is an unsatisfactory situation, not only 
for general reasons but also in view of the rapid experimental progress. 

The canonical way to investigate the ultraviolet divergencies is to calculate 
the corresponding heat kernel coefficients. For the dielectric (nondispersive) ball 
this had been done in Jl| and for the dielectric cylinder in ||, where it had been 
shown, for instance, that the coefficient is zero in dilute order and nonzero 
beyond. In the present note we calculate the relevant heat kernel coefficients for 
a dielectric ball with dispersion. 

Dispersion means a frequency dependent permittivity, e(uj). This is motivated 
by the expectation that for high frequency the permittivity tends to unity so 
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that the ultraviolet modes contribute less to the ground-state energy and the 
divergencies become weaker. 

It is reasonable to assume the asymptotic behavior of e(u) to be 

for uj — > oo. Higher terms do not contribute to the ultraviolet divergencies. Let 
us note that this asymptotic behavior is typical for solid state models, the Drude 
and plasma models for instance. In the latter, the parameter e\ in Eq. (HD has 
the meaning of the plasma frequency squared. 

We remind the reader of some basic formulas. In zeta functional regularization 
the ground-state energy is given by || f§] 

£o(s) = ^EA)- 2s , (2) 

Z 3 

where Xj are the corresponding (discrete) energy eigenvalues. It can be expressed 
in terms of the corresponding zeta function, 

Eo(s) = ^as-\). (3) 

Here \x is an arbitrary parameter with the dimension of a mass. Dropping the so 
called Minkowski space contribution the zeta function can be represented as 

= !^z£ g (2/ + 1} r dk k -2 S d ln 

7T J~[ JO OK 

where fi(ik) is the Jost function of the corresponding scattering problem. A 
detailed explanation of these and related formulas can be found in |], ^, 

The heat kernel coefficients can be obtained from the zeta function by means 

of 

a n = (4nf 2 Res r(s)((s) (5) 
s = \-n 

and the divergent part of the ground-state energy in zeta functional regularization 
is given by (for a massless field) 

ndiv ( - \ ~ a 2 /I 



Here we drop contributions from ax/2 and 03/2, as these coefficients will turn out 
to vanish for the problem considered. 

It is known that in the zeta function regularization there is a smaller number of 
singular contributions to the vacuum energy than in other regularization schemes. 
For example, in the regularization 

£oW = ~I>e-^ (7) 
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with an exponentially damping function the divergent part of the ground-state 
energy is 

Et(5) = T ^(f 2 a 1 + lnSa 2 ), (8) 

again with the Minkowski space contribution already subtracted and dropping 
the terms with ai and as. 

2 2 

The Jost function for the problem at hand is known. For convenience, during 
the calculation we put the radius R of the dielectric ball equal to one, R — 1. 
The dependence on R is recovered by dimensional arguments. With the notation 
v = 1 + 1/2, the Jost function consists of contributions from the TE and the TM 
modes, 

f l {ik)=AZ B {ik)*Z M {ik) (9) 

(/ = 1,2,...) with 

A^ E (ik) = e"5 (Ves'e - se') , (10) 
A™(zjfe) = e-l^s'e-se'V (11) 

Here the abbreviations 

s = 8 t {q) = ffiu(q), (12) 
e = e,(*) = ^K v (k), (13) 

are used where I u (q) and K u (k) are the modified Bessel functions. The prime 
denotes the differentiation with repect to the argument of these functions. The 
arguments of the Bessel functions are related by 

q = V~ek. (14) 

The factors in ( |TUD and flTT|) follow from the normalization condition of the 
regular solution of the scattering problem, for details see the example of a square 
well potential in [0. This is of importance since we consider e depending on k. 
We mention that these representations hold in the presence of arbitrary frequency 
dispersion, as has been noted in (see also ||). 

In order to get the residues according to Eq. (|) it is sufficient to approximate 
the Jost function by its uniform asymptotic expansion for large k and v keeping 



z = - fixed. Using the well known expansions [10], we obtain 



\n Mik) = v(2(r)(^fez)-T)(z))-lne) 

+ In (e^s'e — e~*se'^ + In (e~*s'e — e±se' 
= Do + D te + D T m, (15) 
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where the tilde denotes the Bessel functions without the exponential factors. Now 
we use the expansion 

= ! + % + % + (16) 



as well as the known expression for rj, for example i]'(z) = \/l + z 2 /z, and obtain 

(17) 



n l y/T+J-i l / y/T+^-i e\ (VTT^-i) r 

Jj - -e\ -7: h -r e 2 -: - 



V Z" V 

For D TE and D TM we obtain 

'1 + ez 



l VTT 



+ 



D 



TE 



In 



i 

y2\ 4 



1 + 



1 + C) (1 + B) + 



1 + ez 2 
1 + z 2 



[1 + A)(1 + D) 



(18) 



and 
Dtm 



in <; - 



e 2 



l + C)(l+5) + e 2 



+ ( e 2 - e 2 



l + z 2 7 

(l + jjxi + g) 

4z/((l + 2 2 )(l + e^ 2 ))i 



1 + 



i+a)(i +j d) 

(19) 



where we used the same abbreviations for the Debye polynomials as in [[ij. We 
need them in the first nontrivial order only, A = Ui(tq)/v, B = —Ux{t)/i>, C = 
vx{t q )/v and D = —V\(t)/u with t = + z 2 and t q = 1/vl + ez 2 . Inserting 

now the expansion ( [TBI) of e we obtain finally 



D 



1 ex 



TE 



z/ 3 16 (1 + z 2 fl 2 



+ ... 



and 



D 



1 ex z A + 4z 2 + 4 



TM 



+ 



(20) 



(21) 



z/ 3 16^ 2 (1 + ^ 2 ) 5 / 2 

We have to insert these expansions, Eqs. fllTD , (p0|), (f21|), into the Jost function, 
Eq. (H), and the latter, then, into the zeta function, Eq. (|J). Performing there 
the variable substitution k = uz, the sum over v can be carried out which gives 
Hurwitz zeta functions, 



Ch(s; 



i=i 



(22) 



In summary, we obtain 



+ U2s + 2;-)l i,^ 

ei + 2^ 2 + 2 
~8~z 2 (l + z 2 ) 5 /2 



v/TT^ 2 - l e\ (VTTz^ - 1) 



z A 4 z*y/TT 



Now it is easy to extract the heat kernel coefficients using Eq. The 
rightmost pole is at s — 1/2 resulting from the first Hurwitz zeta function. It 
yields 

a x = -—ex. (24) 
The next pole is at s = — |. It results from the integral in the first line of Eq. 

1 _ 2s d VI + z 2 - 1 -1 1 
^ ^ =27i Sr( " 1_S)r(S+ 2 ) - 

Further contributions result from the pole of the second Hurwitz zeta function. 
Taking all contributions together we obtain 

4^ 2 167T , nr .\ 

a 2 = —e\ + — e 2 , (25) 

where a term linear in e\ cancelled out between TE and TM contributions. To 
Oi, fl2"4|), and a 2 , (^3), the TE and TM modes give equal contributions. 

Formulas ( p4|) and (p5|) are the main result of this paper. In order to draw 
more physical conclusions we restore the dependence on R. From [ex] = R~ 2 , 
[e 2 ] = R~ 4 , [ai] = -R and [a 2 ] = -R -1 we get 

ai = — — eii? 3 , a 2 = — e 2 i? 3 H — — e 2 i? 3 . 

Now the divergent part of the ground-state energy reads in zeta functional regu- 
larization 



s) = -^(\ + 2\n^-2)(ej + Ae 2 )v, (26) 



rndiv/ 

L " 327T 2 Vs 



where V = ^-R 3 is the volume of the ball, and in the regularization Eq. (0), 
using the exponentially damping function, it is 



E ° iv{6) = itA-¥ L+ ( tl+ 4e2 J ln 6 ) v - (27) 

These results show that for a dispersive dielectric ball ultraviolet divergencies are 
present in a fashion similar to that for the nondispersive case. 

We add some remarks. 
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1. In both regularizations the ground-state energy energy is divergent. Be- 
cause of a 2 7^ 0, even after the removal of the diverging contributions an 
arbitrariness (e.g., ln/z in Eq. (|26|)) remains in the finite part. 

The model, Eq. (|l|), chosen for the permittivity reflects the physical as- 
sumption that for high frequencies the dielectric ball becomes transparent. 
From the results, Eqs. (p6|) and (0), in particular from the contribution 
of e 2 , it follows that this is insufficient in order to get a satisfactory phys- 
ical interpretation. Although the dependence of the divergencies on the 
dielectric properties is considerably simpler than in Jl| , we are left with the 
same conclusions as for the nondispersive case that for a dielectric body 
substantial changes in the physical context are necessary [||. 

2. In the sense of renormalization one might try to absorb the divergent contri- 
butions into some classical part. For example, in the bag model, ultraviolet 
divergent contributions proportional to the volume like Eqs. (|2BJ) and (P7|) 
can be put into a redefinition of the bag constant. However, in the present 
case we do not have any classical energy which could be associated with 
the dielectric ball. In addition, we don't have any normalization condition 
which might help to fix the arbitrariness. 

3. As compared with the nondispersive case, [1[], the non vanishing of a 2 is a 
common feature. The only known exception is the vanishing of a 2 in the 
dilute approximation, i.e., to order (e — l) 2 for e — > 1, which allowed the 
ground-state energy to have a physical meaning and ensures that the results 
of different calculations coincide. However, as it was shown in [^] and |rj, 
this is a peculiarity of a ball with sharp boundaries. For a dielectric body 
with non sharp boundaries, i.e., with the permittivity e(r) being a smooth 
function of the radius, a 2 is non zero even in the dilute approximation. 

4. One might hope that the pressure (force per unit surface) is ultraviolet 
finite rather than the vacuum energy itself. For this end one has to divide 
by the surface (47ri? 2 ) of the ball and to take the derivative with respect to 
the radius. As the divergent part of the ground-state energy is proportional 
to the volume of the ball the pressure contains a divergent constant. 

5. One may assume the dielectric ball to be an idealization of a number of 
polarizable atoms. Then a change in the radius leaving the number of atoms 
fixed requires a change in e according to (e — 1) V = const as discussed, for 
instance, in [[11], [12], In this case, by means of Eq. ([]]), ei^V = const 



follows. Due to the presence of ef in (^) and (|27f) , again, a divergent 
contribution is present. 

6. An investigation similar to the present one had been recently carried out 
in |T4| , where the divergent part of the Casimir energy had been calculated 



6 



for the plasma model in zeta functional regularization. This is equivalent 
to calculate the contribution of 6\ to the heat kernel coefficient and is in 
agreement with Eq. (EH). 



7. In the perturbative approach described in |l lj divergencies proportional to 
the surface area were found. It would be very desirable to understand the 
origin of the different predictions in the different schemes used. 
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